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Abatract 

The  use  of  oontioaooa  baaaa  of  representation  other  than  a 
plane -rare  baaia  la  oonaidarad  in  tha  theory  daralopad  originally  by 
Thoaaa  and  Fend.  Tha  baaaa  oonaidarad  hare  are  tha  aata  of  eigen¬ 
function*  of  Haadltonlana  corresponding  to  a  particle  subjected  to  a 
field  of  foroe  which  rarios  imreraely  as  tha  cube  of  the  dlstanoe  to 
aone  fixed  point.  Variation  of  the  strength  of  tha  interaction  rariaa 
tha  baaia  continuously.  Calculations  of  tha  energy  of  the  hydrogen 
aton  ia  carried  out,  with  results  that  are  appreciably  closer  to  the 
quantum  mechanical  result  than  ia  obtained  with  the  original  Thomae- 
Fermi  theory*  luneri eal_resulta  suggest  the  possible  existence  of  a 
statistical  analogue  of  the  KayleiglHtits  rarietion  principle* 


r  » 
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1*  Introduction 

The  statistical  theory  originally  introduced  by  Thomaa"  and 
(o) 

Fermi  *  for  describing  the  behavior  of  msny-electron  systems,  had  its 
origin  in  a  quasi-classical  statistical  mechanical  viewpoint  which  incor¬ 
porated  the  Pauli  exclusion  principle*  Subsequent  investigations  dealing 
with  this  theory  have  provided  it  with  a  better  relation  to  quantum  mechani¬ 
cal  theory  than  it  had  originally*  As  a  result,  several  modifications  of  the 
original  Thomas-Fermi  theory  exist  which  incorporate  features  of  a  dis- 
tinctly  quantum  mechanical  character.  Of  all  of  these,  the  "Inhomogeneity 
correction"  of  von  Weissacker^  has  the  most  significant  corrective  influ¬ 
ence  upon  the  results  obtained  with  the  use  of  the  statistical  theory* 

Although  this  addition  to  the  Thomas-Permi  theory  was  introduced  originally 
(O 

in  an  ad  hoc  manner,  its  presence,  with  slight  modification  from  the 
original,  seems  assured  from  generalised  treatments  of  the  statistical 
method*^ 

However,  in  spite  of  any  improvements  which  have  been  made  on 

the  original  Thomas-Fermi  theory,  there  remains  one  glaring  deficiency  in 

the  theory:  the  density  of  electrons  is  either  singular  in  the  vicinity 

(7) 

of  a  nucleus  or,  on  the  basis  of  somewhat  arbitrary  auxiliary  conditions, 

(8) 

is  set  equal  to  sero*  '  It  is  to  this  aspect  of  the  theory  that  the  pre¬ 
sent  paper  is  directed*  The  use  of  the  momentum  basis,  traditionally 
associated  with  Thomas-Fermi  theory,  is  a  consequence  of  particular  parti¬ 
tioning  of  the  Hamiltonian  operator*  Here  we  address  ourselves  to  the 
problem  of  alternative  partitionings,  preserving  the  spirit  of  the  original 
theory  by  limiting  ourselves  to  those  partitionings  that  lead  to  continuous 


bases 


2.  Form*!  Theory 

In  terms  of  a  general  formulation  of  the  statistical  theory  based 

(9} 

upon  the  use  of  the  density  matrix  of  eany-eleotron  sys  terns ,  the  original 

Thomas-Fermi  density  nay  be  derived  entirely  from  the  aesnaiption  that  the 
kinetic  energy  and  potential  energy  of  a  particle  are  comm  table  observables. 
This  assumption,  in  fact,  is  basically  the  one  cade  by  both  Thome  and  Fersii. 
Such  an  assumption  is  illustrative  of  a  certain  type  which  generates  a  set 
of  quaal-claaelcal  approximations.^10^  For  these  appraxi  mat  ions  the  Hamil¬ 
tonian  operator  of  the  system  is  divided  into  two  parts,  the  choice  of  the 
parts  being  a  matter  of  convenience.  The  assumption  that  the  two  parts 
commute  immediately  produces  an  approximation  to  the  density  matrix  which 
is  the  generalised  analogue  of  the  Thomas-Fermi  theory. 

In  the  interest  of  brevity,  the  formalism  of  (I)  will  be  employed 
here.  The  density  matrix  (representative)  of  a  many-particle  system  then 
may  be  expressed  as^10^ 

■  f*  {«<V2>  *  <«•« 

n-j. 

where  H  is  the  Hamiltonian  operator  of  the  system,  i.41  n(x)  ^  is  any  properly 
symmetrised  (i.e.,  with  respect  to  exchange  of  identical  particles)  complete 
orthonormal  set  of  functions  of  the  correct  boundary  conditions,  x  stands 
for  the  entire  set  of  configurational  and  spin  coordinates  of  the  particles, 
and  Xjj  is  a  parameter  determined  by 

M-Tr  fdxj>M(x,x)  «£  <1'b1«(A„-H)|t|>b>  •  (2*2) 

^  n«l 

The  spectral  operator  0  (Xjj-H)  can  be  expressed  most  simply  in  tens  of  its 
matrix  elements  which  are  evaluated  in  a  basis  of  eigenfunctions  of  H.  Then, 
we  have  for  such  a  set  {  $B(x) }  which  satisfies 


(2.3) 


-it- 


<<p  I  H  I  4?  >  ■  B  S 

NTi 1  -  '  ^n'  n  n,n 

the  matrix  elements 

■  e<VV  *.,n 

with 

0(y)  ■  i»  y  >  o ,  "1 
®(y)  -  o,  y  <  o  .  ' 

In  these  terns  the  "density"  of  the  nany-partlcle  ays  ten  in  its  configuration 
and  spin  space  is  given  by  fH(x,x),  bat  corresponds  to  the  density  asso¬ 
ciated  with  the  M  lowest  energy  eigenstates  of  the  system/^ 

The  procedure  to  be  exploited  here  involves  making  approximations 

to  the  Laplace  transform  of  the  density  matrix  defined  in  Eq.  (2.1).  For 

(12) 

convenience,  we  may  deal  formally  with  the  spectral  operator  and  obtain' 


^e(X-H)^*  (  dX  e"*X  e(A-H),  Re(*)>0 

o 

-*H 

.  5 — :  . 

2 


(2.6) 

(2.7) 


As  a  consequence,  a  representation  for  the  spectral  operator  is  obtainable 
in  terms  of  the  inverse  Laplace  transform,  vis*, 

e(A-H)  -  2I1  e*^  •  *”  »  Re(/ )  >  0  .  (2.8) 

t~U 

Approximations  which  nay  be  made  for  exp  (-JtH)  give  rise  to  approximations 
to  e(X-H). 


In  particular,  the  qaasl-elassical  approximation  to  which  reference 


has  been  made  arises  as  follows.  Let. 


H-Ho+jJiU)  (2.9) 

and  set 

exp  (-* H)  ■  sap  (-  IH^)  exp  (-a H^)  .  (2.10) 

If  m  work  directly  in  terns  of  a  complete  orthonormal  basis  of  properly 
symmetrised  eigenfunctions  of  H^,  say  (  ^°(x) }  ,  which  satisfy 

<< iioi *.,»  <2-u> 

we  obtain  the  quasi-classical  approximation  to  the  density  matrix 

[phk-vl-  , 

IT.  n"1  (2.12) 

1  fa*,  -  <  -  Sl'V>  *  .<**>}  • 

The  preceding  analysis  is  general  and  can  be  taken  over  in  a  single¬ 
particle  approximation  for  a  many-eleotron  system.  In  such  an  approximation, 
is  an  appropriate  single-electron  Hamiltonian  and  is  a  suitable  inter¬ 
action  potential.  The  quantity  M  is  taken  equal  to  the  number  of  single- 
electron  eigenstates  which,  on  account  of  spin  degeneracy,  is  usually  equated 
to  half  the  number  of  electrons  of  the  system. In  the  Thomas -Fermi 
theory  is  identified  with  the  kinetic  energy  operator,  |  <Q  °(x)  ]  is  a 
free  particle  basis  and  H^(x)  is  a  self-consistent  potential  which  is  deter¬ 
mined  from  Poisson's  equation  and  certain  physically  plausible  boundary 
Conditions  upon  the  potential.  Regarding  the  latter,  the  condition  that  the 
potential  has  singularities  corresponding  to  an  attractive  Coulomb  field  at 
any  point  immediately  leads  to  a  singularity  there  in  the  density.  By  con¬ 
trast,  it  will  be  noted,  quantum  mechanical  solutions  yield  no  such  behavior 
for  the  density.  As  a  result,  the  choice  of  a  basis  will  be  examined  with 
the  object  of  alleviating  this  difficulty  in  a  quasi-classical  approximation. 
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3*  The  Basis  Functions 

Because  the  object  of  the  present  paper  is  mainly  constructive 
rather  than  analytical,  no  extensive  analysis  will  be  made  of  the  seat  of 
the  aforementioned  disparity  in  the  density  behavior.  However,  we  m^  take 
note  of  the  fact  that  for  the  systems  generally  of  interest  the  binding 
energy  is  a  finite  quantity.  This  may  be  stated  in  mare  quantum  mechanical 
terms  as  corresponding  to  a  lowest  finite  eigenvalue  of  the  Hamiltonian  of 
the  system.  Mien  measured  relative  to  its  lowest  value,  the  eigenvalues  of 
the  Hamiltonian  are  non-negative  as  is,  thus,  the  (modified)  Hamiltonian 
itself.  However,  in  the  quasi-classieal  approximation  which  leads  to  the 
Thomas -Fermi  theory,  the  Hamiltonian  of  a  single  particle  is  represented  as 
the  sum  of  two  operators!  the  first  is  the  kinetic  energy  and  is  non¬ 
negative  while  the  second  is  the  potential  energy  of  the  electron  and  is 
non-positive  and  can  become  unbounded  for  systems  of  interest.  As  a  conse¬ 
quence,  while  the  Laplace  transform  of  the  spectral  distribution  operator, 

Eq.  (2*7),  is  a  bounded  operator  in  general,  the  approximation  employed  in 
Eq.  (2.10)  does  not  conform  to  this  behavior. 

How,  regardless  of  the  choice  of  a  non-negative  in  Eq.  (2*9), 
it  is  apparent  from  Eq.  (2.12)  that  any  negative  singularity  in  at,  say, ^  ^ 
yields  a  divergent  value  for  the  density.  To  avoid  this  behavior,  H^ 
must  be  chosen  also  to  be  bounded  from  below  (i.e.,  >-*•).  For  systems 
which  involve  only  electrons  and  nuclei,  the  interaction  between  the  par¬ 
ticles  involves  only  Coulomb  terms  which  exhibit  singularities  that  are 
simple  poles.  In  any  approximation  of  a  self-consistent  sort,  we  may  sup¬ 
pose  that  no  singularity  of  greater  order  need  be  introduced.  Hence, 
specialising  to  a  single-particle  approximation,  we  wish  to  examine  the 
possibility  of  finding  ^  and  in  Eq.  (2.9)  such  that  both  are  bounded 
from  below.  That  is,  letting 


•7' 


H  -  T  ♦  V  , 


(3.1) 


vhere  T  Is  the  kinetic  energy  operator  for  an  electron  and  V  its  potential 
energy  function  depending  only  upon  its  configurational  coordinates,  we  seek 
a  potential  function  such  that 


H  *  T  +  V 
— o  —  —o 

Si-2-Io 


(3.2) 


are  both  bounded  from  below.  Since  V  has  no  singularities  other  than  simple 
poles,  it  is  apparent  that  if,  in  the  vicinity  of  the  singularity 


»  «W-— —  ,  ■  >1  (3.3) 

I?  -w 

with  c  a  positive  quantity,  will  be  bounded  from  below.  Similarly,  if  V0 

has  the  behavior  exhibited  in  Eq.  (3*3)  with  s  <  2,  Hq  will  also  be  bounded 
(lO 

from  below.'  It  is  evident  that  a  range  of  choices  is  possible  for  both 
c  and  s  for  the  purpose  at  hand.  Hence,  considerations  of  convenience  may 
be  included  in  any  final  choice  of  V*. 

In  order  to  maintain  some  similarity  to  the  original  Thomas -Fermi 
theory,  it  is  appropriate  to  deal  with  an  Hq  which  admits  only  a  continuum 
of  positive  eigenvalues.  Such  eigenvalues  correspond  to  states  for  which  the 
particle  is  unbound  and  is  reminiscent  of  the  free-particle  basis  of  the 
Thomas  -Fermi  theory.  At  the  sane  time  it  is  desirable  to  have  a  V©  for  which 
the  minimum  value  of  has  as  small  a  magnitude  as  possible.  Because  the 
potential  energy  function  V  is  to  be  kept  general  and  unspecified,  except 
for  its  singularities,  a  condition  leading  to  fulfillment  of  this  desire  is 
impossible  to  formulate  explicitly.  Nevertheless,  a  certain  choice  of  e  and 
s  in  Eq.  (3*3)  will  adequately  fulfill  the  requirements  which  have  been  men¬ 
tioned. 


•8* 


Consider  the  central  field  Hamiltonian 


H 


id  .gf  Sj 


(3.li) 


The  value  of  e  in  Eq*  (3*l»)  n ay  range  from  sero  to  unity.  The  letter  value 
la  the  maximum  which  an  attractive  inverse  square  potential  nay  have  if  the 
ground  state  of  the  ays  tee  is  to  be  unbound*  Any  slight  increase  in  the  value 
of  the  parameters  yields  a  ground  state  energy  value  which  corresponds  to  a 
"fall"  of  the  particle  to  the  origin, i*e*,  a  negative  infinite  value  of 
the  energy*  Thus,  the  systen  corresponding  to  e  ■  1  in  the  Hamiltonian  of 
Eq*  (3*lt)  is  the  limit  ofwnbounded  systems  having  an  attractive  center 
varying  inversely  as  sons  power  of  the  distance  from  that  center.^  With 
c  ■  0,  the  free-partiole  basis  results* 

Each  Schroedinger  equation  corresponding  to  Eq.  (3*U)  generates  a 
complete,  orthogonal  and  continuous  set  of  functions  which  normalise  accord¬ 
ingly.  Because  of  spherical  symmetry  of  H  ,  the  Schroedinger  equation  is 

-*"0 

separable  in  terms  of  spherical  polar  coordinates*  Hence,  if 

S.  ■  hi.  «.5) 


and 


+K4m  •  'fxa(r)  0,1,2,..., 


(3*6) 


where  1^(9,^)  is  a  normalised  spherical  harmonic  of  the  Bulerian  angles 
(17) 

8  and  t  we  obtain  ths  radial  equation 


d2  fggOO  2  dflA  ^  /-2  A  x  1 

*  r  — g -  +  <S  -“jr*  T 


KA 


(r) 


(3.7) 


with 


**■  ft 


(3 .8) 


WLth  transformation 


we 


K  »  Kr 


d2[aft(*) 

- -2 -  ♦ 


(l 


--4->  [r4<r>1 


-  0 


(3.10) 

(3.11) 


which  is  Baa sal's  aquations*  Its  solutions  yield 

^(8)  -  tml/i  Jn  (R)  ,  (3.12) 

JL 

\  m  *  ^(  i  ♦  I)2  -  o/lt 7  (3.13) 

and  Ju(x)  ia  a  Bassal  funotion'  ’  of  tha  first  kind  of  ardor  a. 

For  aach  £  i  0  thera  is  ona  ralua  of  w^  which  is  positive  and 
one  which  ia  negative*  Tha  latter  correspond  to  solutions  which  diverge  at 
the  origin,  while  tha  forwer  do  not*  Moreover,  they  diverge  so  strongly  that 
tha  integral  of  tha  square  of  tha  eigenfunction  ever  any  finite  volume  enclos¬ 
ing  tha  origin  also  diverges*  Henoe,  they  comprise  unacceptable  eigenfunc¬ 
tions  .  Thus,  wa  shall  tske  tha  positive  sign  in  Eg*  (3*13) *  For  o  /  0  tha 
eigenfunction  corresponding  to  £  ■  0  also  diverges  at  the  origin.  However, 
the  integral  of  its  square  over  any  finite  volune  enclosing  it  does  not, 
so  we  retain  it*^^ 

Since  the  Bessel  functions  in  lq*  (3*12)  are  not  quadratieally 
integrable,  wetpaaine  the  integral^ 


■10- 


ahere  I  (x)  la  a  Baanl  function  of  the  first  kind  of  purely  1ml  gin  ary 
u 

(21) 

argument.  From  the  asymptotic  properties  of  the  latter,  the  dominant 

hh 

e 

Wn 

\  lim  5 - ^J-t—  J  .  (3.11*) 

L  .2^0  (**2)*  J 

The  term  in  brackets  is  a  well-known  representation  of  a  delta  function*  As 
a  consequence, 
ee 

(dl^  Ij2  lia  "  1*  *2  0  ;  (3.15) 

o  a*  -qmO 

and  the  original  functions 

+«.  -  -^1  Vb)  1«.-<e,<8>  OM) 

(Sr)*  ’ 

are  properly  normalised  eigenfunctions  of  Sq.  (3*5)  in  the  sense  of  the 
Integration  of  Sq.  (3*16)  which,  together  with  Eq*  (3*7),  preserves  the 
similarity  with  an  integration  over  momenta  as  in  the  Thomas-Fermi  basis. 

It.  The  Quasl-Classlcal  Density  Matrix 

The  quasl-classical  single-particle  approximation  to  the  density 
matrix  now  may  be  expressed  in  terms  of  the  basis  which  has  been  described. 
Making  use  of  Bqs.  (2.12),  (3.2*),  (3.8),  (3.15)  and  (3.16),  we  obtain 


term  is  given  by 


x  (5h)J-*y 

uV  2a  U  KXK 2J 


so  that 


“•  VW  '  irj 


a2  -VO 


*ih 


■11- 


1T  c 

*  BT75 

rl 


-  V(ri» 


(Kri) 

-a - i-  r.  (©. 

i  *»  - 

(Kr^4 


(h.l) 


2  k2  <K2  c 

"BT  ~5 

r2 


Jttp  (KrJ  ?~| 

-»<**»  X  ♦,)}]. 


(&2)! 


With  the  aid  of  the  Addition  Theorem  of  Legendre  polynomials,  '  the  previous 
equation  reduces  to 


0».2) 


where 


°j  m  e(XM  "^ST  75  “  v(rj))  • 

PJ 

The  egression  in  Sq.  (It. 2)  can  be  reduced  further.  However,  we 
shall  only  do  so  for  the  diagonal  elements  of  the  density  matrix,  which 
correspond  to  the  density.  Making  use  of  the  properties  of  the  6-function, 
we  obtain 

oo  Ijf 


Z  (2  0  +  1)  f  dKKJ, 
U*r  tm  o  J 


(Kr) 


(14.3) 


where 


K^j2  ■  (^jj  “  7(r))  ~  cAr2  ^  0  . 


a*.i») 


Apart  from  the  difference  in  the  orders  of  the  Bessel  functions,  the  presence 
of  the  inverse  square  potential  and  a  neglect  of  spin  multiplicity,  Sq. 

(U.3)  la  formally  identical  with  the  corresponding  expression  obtained  in  a 


•12“ 

(23) 

fr «8 -par tide  basis.  Thus  the  terms  corresponding  to  eaoh  Talus  of 

represent  the  contribution  to  the  density  Bade  by  particles  having  a  square 

p 

angular  momentum  of  £  (  £  +  1)  4  • 

The  approximation  to  the  density  given  in  Eq.  (It  .3)  can  be  seen  to 
exhibit  no  divergence  at  the  origin  providing  7(r)  exhibits  no  divergence  as 
strong  as  l/r^.  In  eases  of  atonio  system,  in  which  Coulomb  forces  apply, 
this  condition  is  satisfied*  Because  Kg  vanishes  then  at  som  non-aero  value 
of  t  the  density  vanishes  for  all  smaller  values.  While  this  behavior  is 
clearly  unsatisfactory,  it  possibly  represents  a  better  approximation  to  a 
finite  density  than  is  afforded  by  the  original  Thomas-Permi  theory.  For 
atomic  systems  we  may  anticipate  that  ^g  is  negative.  Thereupon,  an  asymp¬ 
totic  vanishing  of  7(r)  as  r  -►  «*o  indicates  the  existence  of  a  finite 
distance  beyond  which  the  density  vanishes  identically.  Thus,  the  density 
given  by  Sq*  (It .3)  may  be  anticipated  to  have  non-sero  values  in  a  limited 
region  of  space. 

5.  The  Energy 

The  approximation  which  has  been  given  for  the  quasl-clasalcal 
density  matrix  enables  one  to  obtain  an  approximation  for  the  energy.  The 
simplest  way  to  do  so  is  with  the  aid  of  the  relation^ ^ 

M  *M 

<E?g  -  j"XgdM-*gM-J  MdX?  (5.1) 

0  *o 

where  OXg  is  the  sun  of  the  first  M  eigenvalues.  Since  the  quasi- 
claaaical  approximation,  Eq.  (li*2),  involves  ^  g  through  the  spectral 
operator,  we  need  to  examine 


-13 


-57- 7t'»  7«} 

A0  r 

1  9<Vtt-1;  £-7m)- 


When  combined  with  Eqa^  (lj.2)»  (U *3)  and  (5.1),  there  results  the  approxima¬ 
tion 


<*>*  -  /  ^  <2*- pr  *  T<r>}  /" 


ee  f  t» 

f  dr  r  2 
oJ  U-° 


2  <»*♦«!;  /"  *(irt|. 


(5.2) 


The  last  integral  is  the  approximate  expectation  value  of  given  by  Eq. 
(3>i»)  for  the  first  M  eigenvalues  of  H.  Both  integrals  can  be  reduced  con¬ 
siderably  with  the  aid  of  the  Schafheitlin  reduction  formula  for  Bessel 
functional  ^ 

The  present  approximation  gives  no  indication  of  the  reliability 
of  the  energy  values,  there  being  no  extremum  principle  available  to  serve 
in  this  connection.  However,  since  the  eigenvalues  satisfy  the  relation 

<**>,  -  <#  *  {<.>„  -  <«>/  , 


estimate  of  the  precision  of  tbs  approximate  energy  values  can  be  made. 
For  this  purpose,  we  must  have  an  expression  for  the  sum  of  the  squares  of  the 
first  M  eigenvalues.  We  have  ^ 

•  J"  V  «  -  f  XMdA 

0  *0 


(5.3) 


Making  use  of  2q.  (l*. 2)  and  carrying  oat  the  requisite  operation*,  we  obtain 

(  dp  r|  2  Ul*  1){^  ^.»(r)}2  J 

O  IT  0 

o*  r  M  Kg  — . 

'  /arr[t:E  <st*l>^{s  ?  *7(r)'  J 


*  [ 


V  1*  /*  1 

Z  (2Z+  1)  X  f  «JKK5J  2(Kr)  L 

t-0  l»2  J  %  J 


(5.1*) 


With  Eqs.  (5*1)  and  (5»U)  we  are  able  to  estimate  the  precision  of 
the  approximate  energy  values  •  Thus,  for  the  lowest  eigenvalue,  M  ■  1  and 
A  -  -  <i>*  is  a  non-negative  quality  which  vanishes  for  an  eigen¬ 

value.  In  the  present  approximation,  4  Is  a  function  of  c,  the  strength 
of  the  inverse  square  potential.  Hence,  we  may  anticipate  that  some  value 
of  o  may  possibly  render  ^  a  minimum.  For  such  a  value,  the  most  precise 
value  of  <«>x  is  obtained  within  the  framework  of  the  present  theory. 

6.  Computational  Techniques 

The  integration  indicated  in  (1*.3)  may  be  performed  explicitly, 

giving 

«o  _ 

^  «*«♦ «[  y®  -  v<s  • 

i-y,  (6.D 

(l*.l*)  and  u^  by  (3.11*).  This  expression  is  convenient 
for  computational  purposes.  Obvious  alterations  give  density  expressions  for 
the  various  energy  terns  considered  in  the  last  section. 


where  Kg  is  given  by 
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The  definition  of  deneit7  thus  derived  vu  applied  to  the  problem 
of  the  hydrogen  atom,  for  which  7(r)  -  -  efyr.  An  IBM  7090  Bata  Prooesalng 
System  was  programmed  to  evaluate  (6.1),  together  with  Ot.lt)  and  (3.11t),  for 
various  values  of  and  e.  In  all  cases,  summations  over  A  were  extended 

Q 

until  further  contribution  was  less  than  1  part  in  10°  (£~7)*  From  the 
densities  thus  tabulated,  normalisation  and  energy  integrals  were  evaluated, 
using  Slaps on's  Rule.  In  each  case,  the  egression  under  consideration  was 
tabulated  at  a  sufficient  number  of  points  so  that  further  increase  caused  no 
significant  change  in  the  result  of  the  integration.  A  typical  mesh  density 
was  50  points/Bohr  radius.  The  energy  expressions  were  evaluated  and  inte¬ 
grated  in  such  a  way  as  to  allow  the  use  of  some  of  them  as  expectation  values 
of  inverse  powers  of  r. 

Several  sets  of  calculations  were  performed.  The  equation 
M(X  M,c)  »  Mo  was  solved  far  AM  at  various  values  of  c  for  M#  ■  1  and 
Mq  -  5,  corresponding  respectively  to  a  filled  Is  level,  and  filled  Is,  2s, 
and  2p  levels,  and  energies  and  expectation  values  of  powers  of  T  were 
found  for  these  values  of  and  c.  In  addition,  points  on  the  curve 
M(Xm,o)  vs.  X M  for  c  ■  0.5  were  found,  so  as  to  allow  evaluation  of 


(6.2) 


for  Mq  ■  1,  where  A0  is  given  by  M(A0»c  )  "i  ,  the  integration  again  being 
by  Simpson's  Rule.  This  constitutes  a  check  on  the  numerical  procedures,  as 
the  energy  found  in  this  way  must  agree  with  that  found  from  the  integration 
of  the  energy  density  at  A  0  and  c  .  The  number  of  points  on  this  curve,  1*0, 
was  not  sufficient  to  give  a  result  comparable  in  precision  to  the  integration 
of  the  energy  density,  these  calculations  being  used  mainly  to  demonstrate  the 
internal  consistency  of  the  results* 
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7.  numerical  Results  and  Discussion 

Table  IA  presents  the  results  of  calculations  of  the  energy  for 
M  ■  1.  The  most  interesting  feature  of  these  results  is  the  variation  of 
uith  o.  As  can  be  seen,  has  a  definite,  though  broad,  maxlmiwi 

rather  close  to  the  quantum-mechanical  value,  suggesting  the  possible  exis¬ 
tence  of  a  statistical  analogue  of  the  Rayleigh-Ritz  variational  principle. 
Similar  behavior  can  be  seen  in  Table  IB,  which  gives  results  far  M  -  5. 

These  quantities  are  shown  graphically  in  Tig.  la  and  lb.  Table  IA  also 
gives  A  ■  while  a  plot  of  A  vs.  c  is  shown  in  Tig.  2* 

There  exists  a  value  for  c  for  which  A  is  a  minimum,  though  this  does  not 

(57) 

coincide  with  the  "best"  value.'  Also  noteworthy  is  the  great  improvement 
in  the  energy  over  the  value  from  the  Termi-Amaldi  theory,  with  which  these 
results  are  directly  comparable,  since  the  theories  become  identical  for 
c  »  0.  Table  IB  also  presents  expectation  values  of  various  powers  of  r 
together  with  the  Quantum  Mechanical  results  and  Termi-Amaldi  results.  At 
c  -  0.6636,  corresponding  to  the  "best"  value  of  there  is  a  distinct 

improvement  in  all  quantities.  At  e  -  0.95,  approximately  the  value  corres¬ 
ponding  to  minimum  A  ,  the  expectation  values  of  the  positive  powers  of  r 
are  estimated  much  less  well.  Presumably  this  is  caused  by  the  somewhat 
earlier  cutoff  in  density  attendant  upon  the  use  of  c  >  0.  The  removal  of 
the  infinity  in  the  density  it  r  ■  0  is  undoubtedly  responsible  for  the 
improvement  in  the  other  quantities. 

8.  Future  Applications 

So  far,  we  have  constructed  a  generalised  Thomas-Termi  single 
particle  theory.  Application  to  multiparticle  systems  may  be  achieved  by 
invoking  the  idea  of  the  self-consistent  field.  We  may  thus  anticipate 
that  the  improvement  in  energies  found  here  would  be  carried  over  into  a 


self-consistent  field  calculation.  Such  a  aeries  of  calculations  it  in 
progress*  Also,  a  Molecular  system  of  high  symmetry  is  suitable  for  treataent 
in  this  way.  This  work  is  also  in  progress.  Both  will  be  reported  on  in  due 
course* 
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Table  IA 


Sunmaxy  of  Results  for  M  -  1  (in  AtoeLo  Units) 


0 

X 

<A 

A 

<Vr>! 

<r>L 

CM 

0.00 

-.21*038 

-0.7211 

00 

- 

1 .1*1*  2 

OO 

1.51*1 

3.205 

0.50 

-.21*760 

-0.5333 

0.1*1*58 

.1611* 

1.067 

2.1*11* 

1.51*3 

3.117 

0.60 

-.25065 

-0.5289 

0.1,170 

.1373 

1.058 

2.239 

1.528 

3.01,8 

0.6638 

-.25310 

-0.5280 

0.1*01*6 

.1258 

1.056 

2.162 

1.515 

2.992 

0.70 

-.251*70 

-0.5283 

0.3993 

.1202 

1.057 

2.126 

1.506 

2.956 

0.80 

—.2601*0 

-0.5325 

0.3918 

.1082 

1.065 

2.074 

1.1*71* 

2.830 

0.90 

-.2691*0 

-0.51*55 

0.3981 

.1005 

1.091 

2.095 

1.1*21* 

2.61,2 

0.95 

-.27705 

-0.5596 

0.1*123 

.0991 

1.119 

2.161, 

1.382 

2.1*91, 

0.99 

-.28920 

-0.5857 

0.1*1*1*8 

.1018 

1.172 

2.330 

1.319 

2.281 

1.00 

-.30105 

-0.6139 

0.1*81*5 

.1076 

1.228 

2.537 

1.262 

2.096 

Quantun 

-0.5000 

0.2500 

0 

1.000 

2.000 

1.500 

3.000 

at  c  ■  0.50  fron  Integration  of  (6.2)  ■  -0.537. 


Table  IB 


Sunary  of  Results  for  M  ■  5  (in  Atonic  Units) 


0 

X 

<1^ 

0.0 

-0.06221 

-1.2331 

oo 

0.3 

-0.06  2it8 

-1.0711 

0.631*9 

0.1* 

-0.08260 

-1.0606 

0.561*7 

0.5 

-0.08273 

-1.0559 

0.521*2 

0.6 

-0.08287 

-1.051*1 

0.1*933 

0.7 

-0.08303 

-1.0578 

0.1*770 

Quantun 

-1.0000 

0.3125 
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